home *** CD-ROM | disk | FTP | other *** search
/ Aminet 5 / Aminet 5 - March 1995.iso / Aminet / dev / misc / LEDA_gene.lha / LEDA-3.1c-generic / incl / LEDA / matrix.h < prev    next >
Encoding:
C/C++ Source or Header  |  1994-08-05  |  2.3 KB  |  103 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA  3.1c
  4. +
  5. +
  6. +  matrix.h
  7. +
  8. +
  9. +  Copyright (c) 1994  by  Max-Planck-Institut fuer Informatik
  10. +  Im Stadtwald, 6600 Saarbruecken, FRG     
  11. +  All rights reserved.
  12. + 
  13. *******************************************************************************/
  14.  
  15.  
  16. #ifndef LEDA_MATRIX_H
  17. #define LEDA_MATRIX_H
  18.  
  19. //------------------------------------------------------------------------------
  20. //  matrices
  21. //------------------------------------------------------------------------------
  22.  
  23. #include <LEDA/basic.h>
  24. #include <LEDA/vector.h>
  25.  
  26.  
  27. class matrix
  28. {
  29.   vector** v;
  30.   int  d1;
  31.   int  d2;
  32.  
  33.   void     flip_rows(int,int);
  34.   void     check_dimensions(const matrix&) const; 
  35.   double&  elem(int i, int j) const { return v[i]->v[j]; }
  36.   double** triang(const matrix&, int&) const;
  37.     
  38. public:
  39.  
  40.   matrix(int=0, int=0);
  41.   matrix(const matrix&);
  42.   matrix(const vector&);
  43.   matrix(int,int,double**);
  44.  
  45.   matrix& operator=(const matrix&);
  46.  
  47.  ~matrix();
  48.  
  49.   LEDA_MEMORY(matrix)
  50.  
  51.  
  52. int     dim1()  const  {  return d1; }
  53. int     dim2()  const  {  return d2; }
  54.  
  55. vector& row(int) const;
  56. vector  col(int i) const;
  57. matrix  trans() const;
  58.  
  59. matrix  inv()   const;
  60. double  det()   const;
  61.  
  62. matrix solve(const matrix&) const;
  63. vector solve(const vector& b) const { return vector(solve(matrix(b))); }
  64.  
  65. operator vector() const; 
  66.  
  67. int     operator==(const matrix&)    const;
  68. int     operator!=(const matrix& x)  const { return !(*this == x); }
  69.  
  70. vector& operator[](int i)    const { return row(i); }
  71.  
  72. double& operator()(int,int);
  73. double  operator()(int,int) const;
  74.  
  75. matrix operator+(const matrix&);
  76. matrix operator-(const matrix&);
  77. matrix operator-(); // unary
  78.  
  79. matrix& operator-=(const matrix&);
  80. matrix& operator+=(const matrix&);
  81.  
  82. matrix operator*(double);
  83. matrix operator*(const matrix&);
  84. vector operator*(const vector& v) { return vector(*this * matrix(v)); }
  85.  
  86. friend ostream& operator<<(ostream&, const matrix&);
  87. friend istream& operator>>(istream&, matrix&);
  88.  
  89. };
  90.  
  91. inline void Print(const matrix& m, ostream& out=cout) { out << m; }
  92. inline void Read(matrix& m, istream& in=cin)          { in >> m;  }
  93.  
  94. inline int compare(const matrix&, const matrix&) 
  95. { error_handler(1,"compare not defined for type `matrix`"); 
  96.   return 0;
  97.  }
  98.  
  99.  
  100. LEDA_TYPE_PARAMETER(matrix)
  101.  
  102. #endif
  103.